Comprehensive analysis of long-term trends, meteorological influences, and ozone formation sensitivity in the Jakarta Greater Area

Jakarta Greater Area (JGA) has encountered recurrent challenges of air pollution, notably, high ozone levels. We investigate the trends of surface ozone (O3) changes from the air quality monitoring stations and resolve the contribution of meteorological drivers in urban Jakarta (2010–2019) and rural Bogor sites (2017–2019) using stepwise Multi Linear Regression. During 10 years of measurement, 41% of 1-h O3 concentrations exceeded Indonesia’ s national threshold in Jakarta. In Bogor, 0.1% surpassed the threshold during 3 years of available data records. The monthly average of maximum daily 8-h average (MDA8) O3 anomalies exhibited a downward trend at Jakarta sites while increasing at the rural site of Bogor. Meteorological and anthropogenic drivers contribute 30% and 70%, respectively, to the interannual O3 anomalies in Jakarta. Ozone formation sensitivity with satellite demonstrates that a slight decrease in NO2 and an increase in HCHO contributed to declining O3 in Jakarta with 10 years average of HCHO to NO2 ratio (FNR) of 3.7. Conversely, O3 increases in rural areas with a higher FNR of 4.4, likely due to the contribution from the natural emission of O3 precursors and the influence of meteorological factors that magnify the concentration.


Seasonal variation and temporal trends of ozone (O 3 )
Figure 1 depicts the time series of maximum daily 8-h average (MDA8) monthly anomalies and the trends in the 10th, 50th, and 90th percentiles at five sites in urban Jakarta (hereinafter refers to DKI) and one rural site of CBR.All DKI sites generally exhibit a downward trend of surface O 3 in the 50th percentile, as presented in Figure 1.Map of the study area and the time series of the monthly mean anomaly of ozone maximum daily 8-h average (MDA8 O 3 ) during the study period at Jakarta and Bogor sites.Color solid lines denote linear trends in the 90th (blue), 50th (red), and 10th (green) percentiles.Pink shade represents Jakarta Province with sites in central Jakarta (DKI1), north (DKI2), south (DKI3), east (DKI4), and west (DKI5).CBR represents the site in rural Bogor.The shapefile for Jakarta Greater Area boundary map is obtained from Geospatial Information Agency of Indonesia https:// geose rvices.big.go.id/ petar bi/ and the map plot is generated using ArcGIS Pro online 2.8.0.The plot graph is generated using RStudio 1.4.1106.2021.
Meanwhile, MDA8 O 3 anomaly trends in the 90th percentile show different directions from one site to another.Central (DKI1), south (DKI3), and east (DKI4) Jakarta shows a decreasing trend in different order of magnitude.DKI1 and DKI4 experience faster concentration decreases and significant rates of − 10.7 µg m −3 year −1 and − 7.7 µg m −3 year −1 .DKI1 and DKI4 perform at faster-decreasing rates than their lower (50th and 10th) percentiles, indicating reduced extreme ozone episodes during the study period.Despite the reduction rate of O 3 in the 50th and 10th percentiles, DKI2 and DKI5 present insignificant upward trends in the 90th percentile with a rate of 1.9 µg m −3 year −1 and 5.9 µg m −3 year −1 , respectively, suggesting increase number of extreme O 3 concentration.
The MDA8 O 3 (average from all DKI sites) exceedance frequency (> 100 µg m −3 for 8-h threshold) in the dry season showed a sharp rise after 2010, fluctuated until 2018, but underwent a remarkable reduction in 2019, as demonstrated in Supplementary Fig. S1a.The minimum (maximum) exceedance frequency was observed in 2019 (2012).A yearly average of MDA8 O 3 concentration surpassing the 8-h threshold, calculated from daily data, shows a slight reduction from 2011 to 2019 (see Supplementary Fig. S1b).This finding corresponds with the significant downward trend in high O 3 concentration (90th percentile) observed at specific DKI sites, as depicted in Fig. 1 and Table 1.
The MDA8 O 3 increased during the dry period across DKI and CBR sites, peaking in October.Figure 2 displays a map of the average maximum daily 8-h of ozone (AVGMDA8) in urban Jakarta and rural Bogor during the study period for the dry season (April to November).The AVGMDA8 varied and showed consistency with the 1-h O 3 value, with the highest concentration at DKI3, followed by DKI5, DKI2, DKI4, and DKI1 with values of 128, 127, 122, 111, and 99 µg m −3 , respectively.During 10-year period from 2010 to 2019, all DKI sites  www.nature.com/scientificreports/exceeded the National Ambient Air Quality Standard (NAAQS) for 1-h concentration of O 3 (1-h O 3 >150 µg m −3 ) with 18%, 26%, 30%, 37%, and 41% of exceedances, respectively, at DKI1, DKI4, DKI2, DKI3, and DKI5, as presented in Supplementary Fig. S2.The CBR experienced the 1-h O 3 concentration threshold exceedances of 0.1% during 2017 to 2019.The monthly variations in 1-h O 3 at the DKI and CBR sites peaked between May and November with an upward trend.The O 3 concentrations gradually decreased from December to April.During the observation period, DKI5 experienced the highest annual average of 1-h O 3 concentration compared to other sites with a maximum (minimum) of 194 µg m −3 (97 µg m −3 ), followed by DKI2 and DKI3, which were 160 µg m −3 (107 µg m −3 ) and 147 µg m −3 (103 µg m −3 ), respectively.
The overall O 3 level at the CBR site was considered low compared to that in DKI sites and the NAAQS; however, the annual 1-h and MDA8 values increased from 2017 to 2019.MDA8 O 3 anomaly increased in all percentiles by approximately 1.3 µg m −3 year −1 (90th), 4.2 µg m −3 year −1 (50th), and 3.4 µg m −3 year −1 (10th), indicating the enhanced presence of precursors that have critical roles in photochemical O 3 formation in this area.Changing O 3 concentration in the 10th percentile represents changes in the local background of O 3 concentration 24 .

Meteorological impact on the MDA8 O 3 variation
The MLR model yielded the adjusted coefficients of determination (R 2 ) for Jakarta (average from five DKI sites) and CBR areas for all seasons of 0.28 and 0.43, respectively, implying that 28% and 43% of the MDA8 O 3 daily variability were associated with meteorological conditions.The contributions of meteorological and anthropogenic factors to the monthly mean MDA8 O 3 anomalies in Jakarta and CBR areas over the entire period are depicted in Fig. 3a,b, respectively.Meteorological components influenced inter-annual O 3 changes in Jakarta, accounting for 30% of the relative contribution for the entire season of 2013-2019.In general, meteorological components contribution ranged from 0.2 to 98.5% and was not considerably different from that of the dry season (Supplementary Fig. S3a).A slightly lesser meteorological influence shows during the wet season with a 25% relative contribution to O 3 variability.www.nature.com/scientificreports/Pronounced positive contributions from meteorological features such as February 2018 altered MDA8 obs very slightly.Anthropogenic factors contributed 70% to the inter-annual variability of ΔMDA8 obs , averaging 1.5-99.8%.Notably, positive ΔMDA8 obs in 2013, half of 2014, and 2018 were primarily attributed to anthropogenic rather than meteorological.The positive anthropogenic contribution reached more than 90% and was responsible for the increase of monthly mean ΔMDA8 obs during those years.The monthly mean of ΔMDA8 obs fluctuated and increased significantly in 2018 (Fig. 3a).These results indicate that changes in anthropogenic emissions are more favorable for O 3 production, particularly during the dry season.
Downward UV radiation at the surface (UVB) significantly contributed to the O 3 variation in all seasons, although the coefficient was minimal (Table 2).UVB also revealed a positive correlation with the correlation coefficient of r = 0.32, the second-highest value after wind speed (Supplementary Table S2).The UVB pattern was synchronous with maximum temperature (Tmax), particularly during the dry season, as shown in Supplementary Fig. S4c,d.Intense UV radiation is favorable for ozone production.Furthermore, higher temperatures can increase HCHO concentrations from biogenic emissions 25 .More than 70% of the daily observed MDA8 was above the NAAQS and coincided with UVB > 70,000 J m −2 .According to Indonesian Meteorological Office (BMKG), the UV index in Jakarta at 10:00 local time was estimated to be 2-6, indicating a moderate to high risk.A negative anomaly of UVB influenced the O 3 decline in 2013 (− 8.5 µg m −3 ) (Supplementary Fig. S4d).However, the high negative anomaly of UVB in 2016 (− 7759 J m −2 ) was likely canceled out by the positive contribution of Relative Humidity (RH), thus inducing only a small perturbation in O 3 .
UVB and wind speed mean (WS_MEAN) showed higher correlations than other parameters in Jakarta.The wind speed significantly (p = 2 × 10 −16 ) contributed to the stepwise MLR.Wind speed is negatively correlated with MDA8 O 3 .During the dry period, air stagnation frequently occurred in Jakarta, boundary layer height (BLH) was low, calm winds dominated and recorded as approximately 70% of prevailing winds, and the wind speed was < 3 m s −1 leading to pollutant accumulation on the surface 26,27 .There is likely a contribution from the significant negative anomalies of the mean wind speed in 2018 (− 0.2 m s −1 ) to the increased MDA8 met (5.9 µg m −3 ) during the dry period (in Supplementary Fig. S3a).The BLH contributed negatively to Jakarta' s MDA8 O 3 change, indicating that shallow BLH increases O 3 concentration.In addition, the immense positive contribution from emission drivers in 2018 substantially influenced the overall increase in O 3 levels at the DKI sites.High RH is associated with enhanced cloud cover and a greater chance of rainfall, slowing the photochemical reaction of ozone formation 3,4 .
Meteorological and anthropogenic components shared almost equal proportions over 3 years' period of measurement from 2017 to 2019 at CBR site, with average proportions of 52% (ranging from 17 to 91%) and 48% (ranging from 9 to 83%), respectively.However, the number did not change significantly during the dry season.UVB (p = 3 × 10 −6 ), BLH (p < 2 × 10 −16 ), Tmax (p = 0.005), and Surface Pressure (SP) (p = 1.7 ×10 −6 ) significantly contributed to O 3 variation at CBR.The CBR site exhibits complex terrain and thus presents unique weather conditions.High-altitude areas experience lower pressure.With a forced inversion layer and slow wind speeds, vertical mixing and horizontal dispersion in the atmosphere are poor and hindered by mountains, leading to increased concentrations of pollutants 28 .BLH was significant in CBR and had moderate positive correlations.In higher-altitude areas, during the development of a boundary layer, intrusion of O 3 from the upper atmosphere to near the surface is possible, thus contributing to the surface O 3 level 29,30 .This could be a possible reason for a positive correlation between BLH and O 3 concentration.A similar characteristic of BLH and O 3 variations observed at the CBR site is also evident in other locations with comparable geographic features.For instance, Mount Lulin in Taiwan 31 and certain sites in in Beijing-Tianjin-Hebei Province and the Yangtze River Delta, China 4 .

Spatial and temporal analysis of OMI NO 2 and HCHO column density
The NO 2 monthly average in Jakarta ranged from 1.5 to 8.5 × 10 15 molecules cm −2 , with a 10-year average of 3.6 × 10 15 molecules cm −2 (Fig. 5a).This value is lower than those of cities in Japan 13 , China 9,11 , London 32 , and Mexico city 12 for the same study period.The trend in the 50th percentile depicted an insignificant downward trend, with rates of − 0.04 × 10 15 and − 0.01 × 10 15 molecule cm −2 year −1 in Jakarta and CBR, respectively.The tropospheric NO 2 in CBR is lower than in Jakarta, with a 10-year average of 2.9 × 10 15 molecule cm −2 and ranging from 1.1 to Table 2. Coefficient of several variables from the multiple linear regression (MLR) model results for Jakarta (2013-2019) and CBR sites (2017-2019) in the dry, wet, and all seasons.P value coefficient is presented in Supplementary Table S3. ) molecule cm −2 in Jakarta (CBR).10-year average of OMI HCHO column density in urban Jakarta was 12.6 × 10 15 molecules cm −2 with monthly average ranging from 0.9 to 19.4 × 10 15 molecules cm −2 (Figs.4a and 5b).This number is greater than that of urban areas in China (~2 to 18 × 10 15 molecules cm −2 in Beijing and Shanghai), Japan (~1 to 16 × 10 15 molecules cm −2 in Tokyo, Osaka, and Nagoya), and Mexico City (~3 to 9 × 10 15 molecules cm −2 ) 12,33 .
OMI HCHO in CBR area is lower than in Jakarta, with a yearly average of 10.9 × 10 15 molecules cm −2 over the entire study period (Fig. 5b).Biogenic VOC emissions from plants and vegetation are more reactive than anthropogenic VOCs, causing relatively high HCHO columns over cropland areas such as CBR.However, anthropogenic Non-Methane VOC (NMVOC) emissions surpass biogenic VOC emissions 11 in urban cities where anthropogenic activities are massive and intense.Anthropogenic NMVOC emissions in Bogor from on-road motorized vehicles was around 9.7 Gg year −1 in 2016 34 , one-fifth smaller than Jakarta' s emission for the same sector (48.6 Gg year −1 in 2015) 35 .

Seasonal variation of OMI NO 2 and HCHO
Seasonal variations in NO 2 concentrations in Jakarta and CBR exhibited the same pattern, increasing in April and gradually decreasing in October.In Jakarta, the maximum (minimum) NO 2 concentration occurred during the dry (wet) period in July (December) (Fig. 5d,g, and Supplementary Fig. S7).Enhanced column density occurred in northern Jakarta and extended to Bogor city center during the dry period (Fig. 5d).We also examine the monthly variability of NO 2 from in-situ measurements and reveal that the concentration fluctuates throughout the year, and the seasonal pattern is not too obvious, as depicted in Supplementary Fig. S7.Sofyan et al. 36 found no significant differences in fuel consumption in Jakarta between the dry and wet seasons.However, there is an apparent effect of sea breeze on transporting NO 2 from north Jakarta to central and southern Jakarta.On dry www.nature.com/scientificreports/days, converging sea breezes in the southern part of Jakarta cause pollutant accumulation in central Jakarta 36 .
The pollutants then move along the sea breeze front and penetrate the Bogor area.Generally, NO 2 in Jakarta was higher than in CBR for the whole seasons.The maximum NO 2 concentration in CBR occurred in September, and the minimum in December.The reduction in NO 2 was mainly due to meteorological conditions during the wet period.Disregarding exceptional circumstances such as movement restrictions during the COVID-19 pandemic, which were not considered in this study, anthropogenic emissions in Jakarta and CBR were steady throughout the year.
Compared to NO 2 , HCHO was higher during the dry than wet period but fluctuated more (Fig. 5e,h).A slightly higher temperature during the dry period accelerated VOC photochemical oxidation, thereby contributing to high levels of HCHO.Elevated HCHO levels were notable and extended northwest and northeast of Jakarta during the dry period from 2010 to 2019, as shown in Fig. 5e.The HCHO column density in the whole Bogor area varied significantly seasonally.However, it was slightly lower in northwestern Bogor than in central Bogor.

OMI FNR
Figure 5c presents the spatial yearly average FNR value from 2010 to 2019 in JGA.FNR in Jakarta steadily increased to 4.5, consistent with the increase in HCHO column density (Fig. 4a).Seasonal variations in the FNR occurred in Jakarta and CBR (Fig. 5f,i, and Supplementary Fig. S9b).As expected, the FNR decreased in Jakarta' s www.nature.com/scientificreports/dry period but fluctuated slightly in CBR.The reduction in FNR during the dry period was consistent with that shown in Supplementary Fig. S7, where the NO 2 tropospheric column increased.
During the dry period, the maximum 1-h O 3 and MDA8 O 3 showed more days exceeding the 1-h and 8-h NAAQS thresholds than the wet period.It is worth noting that a large positive MDA8 O 3 anomaly in 2018 was mainly driven by a large positive anomaly in anthropogenic activity.The increased MDA8 O 3 is likely due to a response from the increase of OMI NO 2 under a high FNR value of 3.9.However, to confirm this finding, further modeling work should be conducted.

Discussion
O 3 concentration and its formation mechanism are influenced by meteorology and changes in VOC and NO x .The O 3 formation sensitivity could be either VOC-sensitive, NO X -sensitive, or transitional, where an increase of both precursors affects the O 3 production.In a vibrant urban agglomeration such as JGA, characterized by intense emissions and limited and scattered AQMS, investigating the long-term trends of O 3 pollutants and their formation regimes becomes imperative for effective air quality management.Identification of O 3 formation sensitivity aids policymakers in designing appropriate control strategies.The current study evaluates trend and meteorological influences on O 3 concentration variability.Additionally, we leverage space-borne NO 2 and HCHO tropospheric column density to compute FNR as a VOC/NO X ratio proxy to characterize O 3 formation regimes.
Generally, our results show a downward trend of MDA8 O 3 anomalies at all sites except CBR of rural Bogor.Despite the reduction trend, some sites experienced an increasing trend of high O 3 episodes, such as in northern and western Jakarta, and almost steady state in southern Jakarta.This feature leads to high episodes exceeding the national standard for 1-h and 8-h O 3 concentration from 2010 to 2019.As much as 41% of 1-h O 3 concentration exceeded the national threshold for the whole seasons from 2010 to 2019 and 93% for MDA8 O 3 in the 2018 dry season for Jakarta, with annual averages reaching 158.9 µg m −3 .It is worth noting that high positive anomalies of O 3 during 2018 are influenced mainly by anthropogenic drivers rather than the meteorology, as can be seen in Fig. 3a.It is likely due to a response from the increase of OMI NO 2 under increase FNR value of 3.9 (Fig. 4b).However, to confirm this finding, further modeling work should be conducted.
DKI1 is in the city center, while DKI2 is near the busy harbor, power plants, and industrial areas.Both sites are characterized by higher NO X emissions owing to high anthropogenic activities, as reflected by higher NO 2 concentration (see Supplementary Fig. S8).In the vicinity of high NO emissions, O 3 is removed via NO titration.However, its concentration increases further downstream, as in DKI3 and DKI5.DKI3 and DKI5 presented higher exceedance episodes and the least declining O 3 trend.During the dry period, wind direction toward Jakarta mainly originated from the east and southeast, allowing the transport of pollutants from outside Jakarta to the eastern part of Jakarta, where industrial and manufacturing centers are located (around the Bekasi industrial area), contributing to the increase in O 3 concentration at DKI5 (western Jakarta).In July, the HYbrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model' s backward trajectory showed that air mass transport originated from Karawang and passed through Bekasi and South Jakarta (Supplementary Fig. S10).Despite the synoptic scale of easterly and southeasterly winds during the dry season, Jakarta was affected by sea breezes.According to Kitada et al. 37 , the southerly land breeze carried pollutant-rich but ozone-poor air masses to the Jakarta coast at night.In the morning, photochemical reactions produce O 3 in northern Jakarta, and sea breeze circulation develops in the lower layers under weak synoptic winds 37 .The sea breeze, opposite the easterly and southeasterly winds, became dominant owing to the easterly and southeasterly synoptic scales blocked below approximately 1 km by the mountains in the southern part of West Java.The O 3 -rich air mass was transported to southern Jakarta and Bogor as the sea breeze propagated inland up to 60 km from the Jakarta coastline 37,38 .This sea breeze phenomenon affects the convection activity and transport of atmospheric pollutants, which explains the high concentration of DKI3 downwind of Jakarta.
Measurement of OMI NO 2 in Jakarta shows an insignificant decreasing trend and increasing HCHO from time to time.In metropolitan cities in the US, China, Mexico, and Japan, tropospheric NO 2 outnumbers the HCHO.Thus, the O 3 formation is VOC-limited regime, as shown in Table S5.However, that is not the case for Jakarta, where HCHO column leads.As a consequence, the FNR falls in higher value compared to those cities.Study from Lestari et al. 35 shows that the most significant emission load was from CO, NO X , and NMVOC, at approximately 52.9, 143.9, and 48.6 kilotonnes in 2015, respectively.Significant contributors to NO X were the road transportation sector (57%), power plants (24%), industrial sectors (15%), and residential areas (4%) 35,39,40 .Furthermore, by Lestari et al. 35 , the total emissions inventory of NMVOC in Jakarta was comparable to that of NO X in the same sectors and years.Road transport was the major contributor to NMVOC in Jakarta, with a percentage of 96%, followed by industry at 2%, power plants and residential areas at 1%.According to the Ministry of Transportation of the Republic of Indonesia, there was a total transportation movement of 19.63 million movements/day in 2022.Road transport by motorcycles contributed to 88% of the total NMVOC emissions.Jakarta Provincial Agency documented that there were over 15.8 million motorcycles in the city in 2019, approximately a 2.7% increase from 2018.Motorcycles represent the dominant vehicle fleet in Jakarta Province, accounting for approximately 80% of the total vehicles, followed by cars (17%), trucks (3%), and buses (0.2%) 41 .Similar emission types were also found in Vietnam, where motorcycles mainly contributed to the VOC emissions in Hanoi and Ho Chi Minh City.Motorcycles comprise 80% of the total transportation in Vietnam, and the number of vehicles per 1000 population in Vietnam is lower than that in Jakarta 42,43 .Emissions from motorcycles are likely attributable to high HCHO column density in Jakarta.In addition, NMVOCs are abundant in certain areas, such as gas stations 44 and solvent usage industries 45 .
Different features occurred at the rural CBR site.A 3-year measurement period shows that the O 3 concentration is below the government threshold; however, it shows an increasing rate for all percentiles.The growing O 3 concentration value at the rural site indicates the enhanced presence of precursors that have critical  24 .Bogor has a lower population density and more vegetation than Jakarta.The number of vehicles in Bogor is two million, one-thirteenth of those in Jakarta (26 million) 41 .Bogor, as the outskirt of Jakarta, may receive pollutants transported from Jakarta along with the sea breeze; however, this may not contribute substantially to the total emissions of the CBR.Lower concentrations recorded at the CBR were expected owing to lower emission activity, local characteristics, and meteorological disparity.
Air pollution prevention policies have been implemented in Indonesia and locally in JGA for more than a decade.For example, national kerosene-to-LPG conversion for cooking has been proven to reduce emissions 40 .Decreasing particulate matter, sulfate, and nitrate concentrations in Jakarta since 2006 have been used as evidence for the successful implementation of national and local policies, such as Regulation No. 141 concerning the Utilization of Gas Fuel for Public Transportation and Local Government Operational Vehicles 46 .Another successful implementation policy was issuing the Ministry of Environment and Forestry decree number 20/2017 on the usage of EURO IV for gasoline vehicles in late 2018.One current policy is that of the Governor of the DKI Jakarta, Regulation No. 66/2020 on vehicle emission testing, which regulates penalties for vehicles that violate the regulations.
Despite its ongoing implementation of measures, O 3 pollution episodes persist.A more aggressive NO X reduction should be in place in Jakarta.Thorough investigations of the precursors and emission inventory calculation are essential to enhance our understanding of O 3 formation.Moreover, the thresholds for regime classification in Jakarta may require adjustment compared to global FNR thresholds.Further researches, combining modeling and observations, are warranted to identify effective O 3 reduction scenarios for informed policy recommendations.

Summary and conclusion
This study investigated variations and temporal trends in surface O 3 under two different background conditions: urban (Jakarta, 2010-2019) and rural (Bogor, 2017-2019).We assessed the contributions of meteorological and anthropogenic factors and the long-term O 3 formation sensitivity derived from OMI measurements.Despite the downward trend of O 3 anomalies attributed to a slight reduction of NO 2 concentration, Jakarta experiences recurrent elevated O 3 level surpassing both the 1-h and 8-h national thresholds.Generally, anthropogenic drivers are more dominant than meteorological, irrespective of seasonal variations.However, O 3 pollution can be even more severe under changing climate.Further research will focus in modeling simulations to validate the O 3 sensitivity findings presented in this paper and eventually to establish O 3 chemical regime thresholds, which are crucial for devising effective strategies to mitigate O 3 pollution in JGA.

Study sites and air quality data
Jakarta, the center of government, business, and economic activities, contributes the most to the national gross domestic product with over 190 million USD in 2019 47 .We used air quality data from five AQMS in Jakarta (2010-2019) operated by the Environmental Agency of Jakarta Province, one station in Cibeureum Bogor (2017-2019) operated by the Meteorological, Climatological, and Geophysical Agency (BMKG) in collaboration with the National Institute of Environmental Studies (NIES), Japan, as presented in Fig. 1 and Supplementary Table S1.Further local characteristics of Jakarta and Bogor are described in Supplementary Table S6.We used hourly O 3 concentrations and preprocessed the data to eliminate negative values and outliers.The time series of the hourly concentrations were normalized using the z-score following the methodology used by Ren et al. and Sun et al. 24,48 .Points with an absolute z-score greater than 4 (|zt|> 4) were removed from the time series.The daily average was calculated when a station had more than 60% hourly data, and we considered daily data missing if the availability was below this threshold.

Surface O 3 measurement
Jakarta Environmental Agency continuously measures 30-min surface O 3 using HORIBA APOA-370 analyzer.The APOA-370 continuously monitors atmospheric ozone concentrations using a cross flow modulated ultraviolet absorption method or Non-dispersive ultraviolet absorption (NDUV).The NDUV measures ozone O 3 concentration in sample gases using UV light absorption.UV light is emitted into a gas cell where the sample gas flows, and O 3 absorbs UV light proportionally to its concentration.The absorbed light is detected by a photodiode, generating electrical signals.NDUV employs a cross-modulation method with a solenoid valve unit to switch between sample and reference gases, allowing the detector to distinguish AC and DC signals during processing.AC signals represent gas concentration, while DC signals compensate for UV light source aging, ensuring analyzer stability.CBR site in collaboration with NIES measure O 3 concentration continuously using Kimoto OA-787 49 .

Meteorological data
Meteorological data were used to construct the stepwise MLR model.Parameters such as Tmax, rainfall (RRR), WS_mean, and RH of Jakarta and Bogor were obtained from the BMKG.Other meteorological parameters, such as the BLH, SP, and UVB, were downloaded from the European Center for Medium-Range Weather Forecasts Reanalysis (ECMWF ERA5) product on single levels with a spatial resolution of 0.25° × 0.25° and a temporal resolution of 1 h for the corresponding grid, similar to the meteorological stations and AQMS (https:// cds.clima te.coper nicus.eu/ cdsap p#;/ datas et/).

Trend analysis of ozone metrics
To identify O 3 pollution, we use an indicator from the Indonesian Government Regulation No. 22 Year 2021 on the implementation of environmental protection and management for 1-h O 3 concentration (measured from 11:00 to 14:00 local time) and O 3 metric following the Tropospheric Ozone Assessment Report (TOAR), such as the maximum daily 8-h average (MDA8) 51 .This metric is intended to evaluate the impact of O 3 on human health, vegetation, model comparisons, and the characterization of O 3 in the free troposphere 51 .MDA8 O 3 was calculated from 24 running 8-h averages.The 8-h running mean for a particular hour is computed based on the concentration for that hour plus the following 7 h.The daily maximum 8-h concentration for a given calendar day was the highest of the 8-h average concentrations computed for 8-h periods starting from that day.The 8-h running mean was considered valid when at least 5 h of data were available (60% completeness was required).
The average MDA8 (AVGMDA8) is the mean value of daily MDA8 (calculated based on more than 60% of the daily MDA8).Quantile Regression (QR) estimates the trend in MDA8 O 3 concentration.QR is a well-suited technique for detecting heterogeneous distributional changes, which is often the case for free tropospheric and surface ozone as they typically show diverse percentile trends 52 .Trend uncertainty is estimated by doing a moving block bootstrap algorithm to take autocorrelation into account in the trend uncertainty.Trend analysis considers the seasonality of MDA8 O 3 time series by calculating the MDA8 O 3 anomaly.MDA8 monthly anomalies are more accurate than the monthly average in case of missing data 10,48,53 .Monthly anomalies were calculated by subtracting the individual monthly mean from the monthly mean for the same month for the entire study period.Meanwhile, zero anomaly refers to a situation where the monthly MDA8 O 3 is equal to the long-term average value or no deviation or departure from the average conditions.QR analysis is computed using Rstudio with a quantreg package.

Stepwise multiple linear regression (MLR) model
Statistical models were used to infer the contributions of meteorological and anthropogenic factors to the monthly ozone variability.Researchers commonly use the MLR to quantify the effects of meteorological factors on O 3 4,6,9,10,24 .The MLR uses the following equation to reflect the relationship between a quantitative dependent variable and two or more independent variables: The stepwise MLR estimated MDA8 O 3 from meteorological factors (MDA8 met ); where y is the daily observed MDA8 O 3 (MDA8 obs ) in season s (dry, wet, and all seasons) and city c (Jakarta and Bogor); Met k is the selected daily meteorological parameter; β s,c is the regression coefficient in season s and city c; β o is the intercept; and ɛ is the residual.MDA8 obs used to construct the model were averaged from the MDA8 obs at five sites in Jakarta and one for rural Bogor.
A series of steps was conducted to determine the relative meteorology contribution to O 3 variation.First, correlation coefficients were established between MDA8 obs and all meteorological parameters (see Supplementary Table S2).Second, the Variance Inflation Factor (VIF) was calculated to avoid multicollinearity between variables.Variables with a VIF > 10 indicated multicollinearity removal, and those with a 10 were retained.In this study, all the predictors had a VIF < 10 (approximately 1-3) and were retained for the next selection step.Third, stepwise regression was conducted to obtain the best model fit by adding or removing predictors based on Akaike Information Criterion (AIC) statistics.In this step, optimal meteorological parameters were selected and used to establish an MLR model.The optimal meteorological parameters, calculated intercept (β o ), regression coefficient (β s,c ), and adjusted R-squared for Jakarta (2013-2019) and Bogor (2017-2019) in the dry, wet, and all seasons are listed in Table 2, and Supplementary Table S3, S4.The MDA8 O 3 resulting from the stepwise MLR model shows the contribution of the meteorological component (MDA8 met ).The difference between MDA8 obs and MDA8 met is called the residual and is considered an anthropogenic driver attributed to the ozone variability (MDA8 ant ) 9,10,29 .
The total change in meteorological contribution (ΔMDA8 met ) to O 3 concentration can be calculated as follows: where ΔMet k is the change in the k-th meteorological variable.The change in MDA8 ant (ΔMDA8 ant ) was the difference between the change in observed MDA8 (ΔMDA8 obs ) and ΔMDA8 met .Thus, the relative contribution of meteorology to the ozone concentration could be calculated as the ratio of ΔMDA8 met to the total meteorological and anthropogenic factors (ΔMDA8 met + ΔMDA8 ant ).

Ozone formation sensitivity (OFS)
A common approach involves using satellite measurements to determine observed VOC-to-NO X ratios (FNR).FNRs have been widely used to study OFS.The tropospheric NO 2 column reflects NO X emissions at the surface because most NO 2 satellite measurements are within the planetary boundary layer and quickly transform or are removed from the atmosphere (short lifetime) 14,18 .HCHO is an intermediate product of all VOCs, and can be used as a proxy for VOC reactivity 54 .HCHO is a short-lived molecule that is usually found near emission sources.
The FNR threshold for separating the VOC-sensitive regime from the NO x -sensitive regime was derived using various methods for urban areas in the U.S. and China.The FNR threshold can vary by location owing to the influence of precursor emissions, meteorology, and geography 17,55 .This study compared FNR values in Jakarta and Bogor with those reported in the previous studies.The flowchart and methodology used in this study are shown in Supplementary Fig. S11.

Figure 2 .
Figure 2. Map of the average of Maximum Daily 8-h ozone (AVGMDA8) in urban (Jakarta) and rural (Bogor) during the dry season (April to November) over the entire study period.(ArcGIS Online Pro 2.8.0.2021).The shapefile for Jakarta Greater Area boundary map is obtained from Geospatial Information Agency of Indonesia https:// geose rvices.big.go.id/ petar bi/ and the map plot is generated using ArcGIS Pro online 2.8.0.

Figure 5 .
Figure 5. Map of the average (a) OMI NO 2 tropospheric column, (b) OMI HCHO vertical column, and (c) FNR value from 2010-2019 for all seasons; (d), (e), and (f) for dry seasons; and (g), (h), and (i) for the wet season.Transparent areas indicate an FNR value above 5.The shapefile for Jakarta Greater Area boundary map is obtained from Geospatial Information Agency of Indonesia https:// geose rvices.big.go.id/ petar bi/ and the map plot is generated using ArcGIS Pro online 2.8.0.